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Abstract 

Non-relativistic reduction of the S-matrix for the quasifree electron 
scattering process A ( e, e'p ) A — 1 is studied in order to understand 
the source of differences between non-relativistic and relativistic mod- 
els. We perform an effective Pauli reduction on the relativistic ex- 
pression for the S-matrix in the one-photon exchange approximation. 
The reduction is applied to the nucleon current only; the electrons are 
treated fully relativistically. An expansion of the amplitude results in 
a power series in the nuclear potentials. The series is found to con- 
verge rapidly only if the nuclear potentials are included in the nuclear 
current operator. The results can be cast in a form which reproduces 
the non-relativistic amplitudes in the limit that the potentials are re- 
moved from the nuclear current operator. Large differences can be 
found between calculations which do and do not include the nuclear 
potentials in the different orders of the nuclear current operator. In 
the high missing momentum region we find that the non-relativistic 
calculations with potentials included in the nuclear current up to sec- 
ond order give results which are close to those of the fully relativistic 
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calculation. This behavior is an indication of the importance of the 
medium modifications of the nuclear currents in this model, which are 
naturally built into the relativistic treatment of the reaction. 
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1 Introduction 



Electromagnetic probes provide invaluable information about nuclear struc- 
ture. The quasifree process ( e, e'p ) has been used extensively to study proton 
hole states and to determine single particle spectroscopic factors [|l| . This re- 
action is advantageous because the electromagnetic interaction is known; the 
relative weakness of the reaction permits the probe to interact almost uni- 
formly through the entire nucleus, and the first order of perturbation theory 
should provide an adequate description of the process. Coincidence mea- 
surements of the ( e, e'p ) reaction can provide detailed information about 
the single particle structure of the nucleus over a wide range of momentum 
transfer. 

The ( e, e'p ) reaction has been widely studied both non-relativistically 
[0, 3 and relativistically @, 0, and there are some discrepancies between 
the results of these investigations. Both analyses begin with a lagrangian 
which allows for the interaction of the photon with both electrons and pro- 
tons. Non-relativistic analyses involve the reduction of the free electron- 
proton interaction to a form involving two-component spinors for the nucleon. 
This results in a hamiltonian which is expanded in powers of 1/M where M 
is the nucleon mass The resulting interaction hamiltonian is sandwiched 
between Schrodinger wave functions describing the nucleons in order to form 
the nuclear current. Relativistic analyses are based on the Feynman diagram 
for one-photon exchange between the projectile electron and a proton which 
is imbedded in the nucleus. The electrons and nucleons are all described rela- 
tivistically as spin 1/2 objects via the Dirac equation containing appropriate 
potentials [|, |, |]. A long-standing problem in quasifree electron scattering 
has been that the spectroscopic factors extracted from non-relativistic anal- 
yses are smaller than expected from shell model calculations. Spectroscopic 
factors which are found on the basis of the relativistic approach are generally 
larger than those found via the non-relativistic approach 

Several groups have attempted to understand the underlying differences 
between these two approaches. This mainly involved looking at the sensitiv- 
ity of quasifree electron scattering calculations to different optical potentials 
and renormalizations of the continuum wave function 0, ||, . This concen- 
tration on optical potentials was largely a result of the improvement in the 
description of proton elastic scattering observables via Dirac phenomenology 
over the standard non-relativistic optical model description. Boffi et al. 
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have multiplied the non-relativistic continuum wave function by a potential- 

1 /2 

dependent factor {1 + [S (r) - V (r)] / {E + M)}'', where S (r) is the Dirac 
scalar potential and V (r) is the vector potential. This modification essen- 
tially changes the two-component Schrodinger wave function into the upper 
component of the Dirac wave function, while no other change is made in the 
non-relativistic calculation. They find that extracted occupation probabili- 
ties are larger than those obtained from the unmodified non-relativistic anal- 
ysis. The analysis of Udias et al. replaces the non-relativistic bound state 
wave function with the upper component of a Dirac wave function, and the 
non-relativistic continuum wave function is modified by factors of the same 
shape as the factor used by Boffi et al. The continuum wave function in this 
case is generated from Schrodinger-equivalent potentials [|r^. The nuclear 
current operators are obtained in the standard way by expansion to order 
1/M'^. Their "non-relativistic" calculations then involve non-relativistic nu- 
clear current operators surrounded by the upper components of Dirac wave 
functions. With these choices little difference is found between the relativis- 
tic and "non-relativistic" calculations. Their conclusion is that differences in 
observed cross sections are not due to non-relativistic reduction, but to the 
choice of optical potential. Jin and Onley have presented a model which 
can take either relativistic or non-relativistic optical potentials while keeping 
other aspects of the calculation the same. They find that different optical 
potentials can change the results by as much as 14%. 

These results clearly demonstrate the variability due to final state interac- 
tions, however, the issue is clouded by the occasional use of upper components 
of Dirac wave functions in a non-relativistic calculation. We believe that the 
essential difference between relativistic and non-relativistic approaches are 
not just in the changes in the optical potentials; these are usually phenomeno- 
logical and equivalent potentials can always be found. Rather the essential 
difference is in the appearance of the nuclear potentials in the nuclear cur- 
rent operators when the relativistic amplitude is reduced to a non-relativistic 
form. Such medium effects on the nuclear currents are absent in the standard 
non-relativistic calculations. 

In this paper we study the differences between the relativistic and non- 
relativistic approaches in calculating the amplitude for the ( e, e'p ) reaction. 
We do this through an effective Pauli reduction of the relativistic transition 
amplitude [|l^. An expansion of the amplitude in powers of {E + M)~^ 
allows us to recover a non-relativistic limit, which matches the standard 
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non-relativistic calculations, with the difference that optical potentials used 
to generate the distorted waves are exactly equivalent to those used in the 
relativistic calculations. The main difference, as mentioned above is that the 
nuclear currents are potential dependent. We compare the two approaches 
and thus explain why they can still give different values for the extracted 
spectroscopic factors, even when equivalent optical potentials are used. 

We introduce the relativistic amplitude for quasifree electron scattering 
in section 0. Section |] outlines the Pauli reduction of the amplitude and 
some of its relevant features are discussed in section The non-relativistic 
limit is discussed in section ^ In section ^ we compare our non-relativistic 
calculations with and without nuclear potentials in the nuclear current op- 
erators, to the results of the fully relativistic calculations. Our conclusions 
are given in section |^. 



2 Relativistic Amplitude 



We consider the one photon exchange model for the ( e, e'p ) process [|], 
in which a photon is exchanged between the incident electron and a target 
proton. The struck proton is detected in coincidence with the final state elec- 
tron. In this paper we are interested in the differences between the relativistic 
and non-relativistic treatment of the hadronic part the ( e, e'p ) reaction. In 
the course of this discussion we do not include the Coulomb interaction in 
the leptonic part of the S-matrix since this will only be important for heavy 
nuclei 0, We will therefore not discuss any nuclei heavier than '^'^Zr in 
this work. 

The relativistic expression for the S-matrix describing the quasifree elec- 
tron scattering process ( e, e'p ) in the distorted wave Born approximation 
(DWBA) is H 



(27r)iV2 
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where Sj^j^{Jb) is the spectroscopic factor and and are electron and 
nuclear currents respectively. M is the nucleon mass and Ec is the energy of 
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the outgoing proton. The electron current is given by 

J^{y)=^e,{y)i'i'eAy). (2) 

where ipei and ipef are the initial and final Dirac spinors for electrons. The 
electron wave functions are taken to be free Dirac spinors and the integration 
at the electron vertex can then be done analytically. This also allows the 
momentum integration for the photon propagator to be done, leaving one 
four-dimensional integration at the nucleon vertex. The nuclear current is 
similarly given by 

Jn i.^) = i^Nf (x) Jn^n, (x) , (3) 
where the nuclear current operator j'^^ is the choice cc2 discussed by de Forest 



F.(«^)y' + ^.'V. (4) 

Fi (g^) and F2 {q^) are nucleon form factors and are functions of the four- 
momentum squared of the exchanged photon, which couples to the nucleons. 
We have = k^—k^j where fcf and fc^ are the momenta of the initial and final 
electrons respectively. The matrix a^'^ is formed from the Dirac 7-matrices 
in the standard way as 



= 1 (7'^7^ - 7'^7^) . (5) 

The integration over d^q in Eq. (|I|) is associated with the propagator of the 
exchanged photon. 

The electrons are described by positive energy Dirac spinors, and the in- 
tegration over coordinates at the electron vertex, in the S-matrix of equation 
in Eq. (|^), can be done analytically to yield a Dirac 5- function giving energy 
and momentum conservation at the vertex so q = kf — kj. The resulting 
(5-function fixes the momentum of the intermediate photon so the integration 
over that momentum is done trivially. The integration over the time coordi- 
nate at the nucleon vertex can then be done to yield a 5-function providing 
overall energy conservation. The S-matrix can then be cast in the form 



M ml 
EcEfEi 



5 (E, + Ej-E,-Ei 



Y: {Jf,JB;Mf,MB\J^,M,) [Sj^j^{Jb)] Zaf^f, (6) 



X 

JbMb 



6 



where Z^^l^ is a function of the initial and final spin projections, momenta, 
etc. Specifically we have at this point 



Z'u^lf = J d'x {k„ x) T^^j^M, (x) exp {tq ■ x) 
where the 4x4 matrix operating on the nucleon spinors is 



7o 



(7) 



(8) 



The four-vector which comes from the electron vertex is 
'Ef + meEi + nie^^/"^ 



2mp 



x{l/2,Uf\ 



2mf, 



Ef + TUe 



7o7 



1 

cr ■ kj 

Ei + rrie 



|l/2,z/i). 



(9) 



and this depends on the energies and momenta of the initial and final elec- 
trons as well as their spin projections. Useful details of the Fock space 
calculations and the expansion of the Dirac wave function in partial waves 
can be found in the paper by Johansson and Sherif [|l^ . When the appropri- 
ate factors of h and c are included, the relativistic expression for the triple 
differential cross section is related to by 



a 



a 



dVt^dVtpdEp 



(27r)^ he _ 



{q^cf PiC 
2Jb + 1 ' "f"' ' ^ 



(10) 



where z/j and Uf are the spin projections of the incoming and outgoing elec- 
trons respectively, while Mb and fi are the spin projections of the bound and 
continuum protons. 

The cross section for quasifree electron scattering in the plane wave im- 
pulse approximation can be written in a factorized form as the product of 
three parts [|l^, a kinematic factor, the cross section for the elementary 
process e + p ^ e' + p', which is evaluated off-shell, and finally a function of 
the energy and momentum of the nucleon inside the nucleus referred to as 
the spectral function. 
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In the following we will discuss results of calculations of the spectral 
function, proton polarization and an asymmetry parameter. The spectral 
function is obtained from the cross section given above by dividing by a 
kinematic factor and the cross section for the elementary process for e + p — > 



e + p. We write [[13 



S (pj = dn.dn^dE. ^ ^^^^ 



EcPc 



da 

dQ,e 



free 

where Ec and pc are the energy and momentum of the final state proton, 
and Pm is the missing momentum, i.e. the momentum of the bound nucleon 
in the initial state. The free cross section is calculated using the nucleon 
current operator of Eq. (^, and is evaluated using the kinematics of the 
quasifree process, i.e. off-shell. Note that the experimental data are divided 
by the elementary cross section ccl of de Forest , while we use the nuclear 
current operator which leads to his cross section cc2, throughout this work. 
We are not concerned with detailed comparison with experimental data in 
this work, so we retain a consistent approach by using the same form for 
the current operator in the calculation of the quasifree S-matrix and the 
elementary process. 

The polarization of the final state proton is given by 



Tm V zV2Ms \v-1/2Mb 



2 



(12) 



We also define an asymmetry parameter in the missing momentum which is 
calculated from the differential cross sections of Eq. (|TU|) as 

. _ (fa {Pm > 0) - d^a {pm < 0) , . 

^ d^a > 0) + d^a (p„ < 0) • ^'^^ 

This asymmetry is similar to the parameter defined by Bianconi, Boffi and 
Kharzeev 



3 Pauli Reduction 

We now perform the effective Pauli reduction ||rT| on the hadronic part of the 
amplitude only; the electrons are treated relativistically throughout. Con- 
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sider the nuclear current of Eq. (j^) above. Using Eq. (^) one can write: 

iKF2{q'^ 



2M 



(14) 



The Dirac spinors may be written in terms of their upper components as 
1 



N 



(T ■ p 



M + E + S(r)-V(r) 



(15) 



where S (r) and V (r) are the scalar and vector potentials, respectively, for 
either the bound or final state nucleons. The energy of the nucleon is E, and 
the associated momentum operator is p. The upper component of the Dirac 
spinor u is related to a Schrodinger-like wave function "^sch by ||16 



u 



D^^Sch where D (r) = 



_ E + M + Sir)-V{r) 



E + M 



(16) 



Note that the two-component wave function "^sch is a solution of the Schro- 
dinger equation used in the ordinary non-relativistic calculations, i.e. con- 
taining central and spin-orbit potentials, but these potentials are "equivalent 
potentials" meaning that they are functions of the vector and scalar poten- 
tials of the original Dirac equation, as well as containing explicit energy 
dependence. For the continuum nucleon the Dirac potentials result in an 
improved description of nucleon-nucleus elastic scattering data |T^, while 
for the bound state the Dirac potentials offer a slightly better description of 
the spin-orbit splitting than those used in earlier non-relativistic calculations 



The relativistic nuclear current of Eq. ( |T4D can, with the help of Eq. (|1 
be written in the form 



X 7" 



cr ■ p 



X 



M + E, + S, ir) - V, ir) 



. ' 2M ^ 

1 

cr ■ p 

M + E, + SJr)-VJr) 



(17) 



D'J' (r) ^ *5ch,6. 
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We now perform an expansion of the object between the braces of Eq. (IT 



The usual representation of the Dirac 7-matrices is used |jT^ to write the 4x4 
operator in terms of 2 x 2 Pauh matrices. The matrix multiphcations are 
performed and a 2 x 2 operator results. The radial function D^^"^ (r) from Eq. 
(|16D and the factor [E + M + S [r) — V coming from both the bound 

and continuum wave functions are then expanded in powers of {E + M)~^ . 
This procedure leads to a sum of reduced nuclear current operators for each 
of the contributing orders: 



J" ix) 



u(0) , .u(l) , .u(2) , 



Sch,b [X) 



(18) 



The reduced current operators can be written in terms of time-like and space- 
like components as 
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and the labels b and c refer to the bound and continuum states, respectively. 
Using the current operators from Eq. 
the S-matrix (|^), we find that for the 
order in {E + M)~^ reduces to 



T9|) up to first order in {E + M) , in 
e, e'p ) reaction the S-matrix to first 



S 



(1) 

fi 



-le 
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Note that we have written this equation in a form in which the integrations 
over the electron coordinate and the intermediate photon momentum have 
not been done. The expansion method does not depend on the plane wave 
approximation for the electrons, and electron distortions could be included 
if desired. The S-matrix to second order in {E + M)~^ is similarly found by 
including the second-order nuclear current as well. 

Note the dependence of the nuclear current operators on the Dirac vector 
and scalar potentials ( through the functions Qi ). This dependence appears 
in all orders of the reduction scheme. Thus as we go to a description in terms 
of the Schrodinger-like wave function for the nucleon, the currents undergo a 
medium modification affected via the nuclear potential. This point is central 
to the present work. We shall discuss the traditional non-relativistic limit of 
the amplitude in section ^; but will concentrate in the following section on 
clarifying the role of the nuclear potentials in the convergence properties of 
the Pauli expansion of the S-matrix. 
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4 Convergence of the Expansion 



In this section we discuss the convergence of the expansion obtained above 
to the fully relativistic calculation. In these convergence calculations, all of 
the factors in the expansion are the original relativistic ones. This is not 
yet equivalent to a standard non-relativistic calculationl The non-relativistic 
calculations are discussed below. 

The calculations of the relativistic S-matrix requires knowledge of the 
Dirac wave functions for the bound and continuum states. For the bound 



state Hartree bound state wave functions are used 181. The continuum 



wave functions for the knocked out proton are obtained using the energy 



and A dependent optical potential of Cooper et al. We restrict our 

discussion to the case of parallel kinematics [|l|. In the diagrams referred to 
in the following discussion the curves are labelled according to their order 
m {E + M)~' for the expansion calculations, and whether or not the Dirac 
potentials are included in the nuclear current operators: dotted curve - first 
order in [E + M)~^ without Dirac potentials; dashed curve - first order 
with potentials; dot-dashed curve - second order without potentials; dot-dot- 
dashed curve - second order with potentials; solid curve - fully relativistic 
calculation. In doing these comparisons we are attempting to clarify the 
convergence of the expansion and the role of the nuclear potentials ( as they 
appear in the nuclear currents ) in the rate of convergence of this expansion. 

Figure 1 shows observables as a function of missing momentum for the 
reaction ^^0{ e, e'p Y^N, 1(a) is the spectral function while 1(b) is the proton 
polarization. The ground state of the residual nucleus, ^^A^, is assumed to be 
a Ipi hole. The energy of the incident electron is 456 MeV, and the kinetic 
energy of the detected proton is fixed at 90 MeV with parallel kinematics. 
The relativistic calculations of the spectral function are fitted to the peak 



of the data |2^; the resulting "spectroscopic factor" is then used in all the 



other calculations for that particular state. ( We adopt this simple fitting 
procedure because our main concern here is comparison between the different 
calculations, rather than a judicious determination of the spectroscopic fac- 
tors. ) Note that the calculations with potentials included converge rapidly 
toward the fully relativistic results in this case, with the curve for the sec- 
ond order calculations being very close the relativistic results over the range 
of momentum transfers shown. In calculating the spectral function, the in- 
clusion of the potential in the first order interaction terms bring the results 



12 



closer to the fully relativistic calculation than the second order without po- 
tentials. It must be stressed that the inclusion of potentials in the interaction 
brings the results close to the fully relativistic results, while the calculations 
without potentials are quite far from the relativistic results and do not show 
a strong indication for convergence to the relativistic result. We have also 
done similar calculations for the same target but leaving the residual nucleus 
in an excited state, as well as using different targets, namely ^^Ca and ^^Zr, 
with the residual nucleus left in both ground and excited states. These cal- 
culations show the same behavior as the calculations shown in Fig. l.The 
above results are of course expected on simple mathematical grounds. The 
essential point however, is to shed light on the role of the appearance of the 
potentials in the nuclear currents. We have seen no evidence that expansions 
that are based on free vertices ( i.e. no nuclear potentials ), will converge 
to the fully relativistic results, even if calculations are done to higher order 
in the inverse of the nucleon mass [0. This will have implications for the 
comparisons with the standard non-relativistic calculations, which we discuss 
next. 

5 The Non-Relativistic Limit 

The expansion of the S-matrix in powers of {E + M)~^ discussed above does 
not quite yield the amplitude used in standard non-relativistic calculations. 
Some care must be taken at this point in the discussion to differentiate be- 
tween the correct non-relativistic limit, and the standard operator used in 
non-relativistic calculations. There are three things that must be done in 
order to obtain the proper non-relativistic limit from the relativistic ampli- 
tude: 

i) The bound state wave function must be normalized to unity. In 
the expansion obtained above, the Dirac bound state wave func- 
tion is normalized to unity and the related Schrodinger-equivalent 
wave function is not. In the non-relativistic calculations it is the 
Schrodinger-equivalent wave function that must be normalized. 

ii) The continuum wave function must be normalized correctly. 
The factors arising from the Dirac field and the normalization of 
the Dirac wave function result in a factor of {E + M) /2E being 
set equal to one to obtain the non-relativistic expression for the 
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cross section, ( this is equivalent to multiplying the right-hand- 
side of Eq. (0) by the inverse of this factor ) . 
iii) Finally, to obtain non-relativistic expressions for the nuclear 
current operators from the relativistic expressions of Eq. (p!9D, 
the nucleon energies ( both continuum and bound ) are set equal 
to the nucleon mass, i.e. E —>■ M. 

It is important to note that these changes still have not yielded the stan- 
dard non-relativistic amplitudes because the nuclear current operators at this 
stage contain the Dirac potentials explicitly. This is an essential difference 
between the relativistic and non-relativistic approaches, and the presence 
of these potentials can lead to large differences in the observables obtained 
via relativistic and non-relativistic approaches. In order to obtain the usual 
non-relativistic expression, the Dirac potentials must be removed from the 
nuclear current operators. When this is done, the non-relativistic equiva- 
lent of the S-matrix of Eq. yields the usual first order non-relativistic 
transition amplitude used by many authors p|, |[. When terms to second 
order are included in the non-relativistic S-matrix, and in the limit of no nu- 
clear potentials, there are some differences between our expression and the 
usual non-relativistic second order S-matrix, which is obtained via a Foldy- 
Wouthuysen transformation of the interaction between electrons and free nu- 



cleons P] . Fearing, Poulis and Scherer have compared Foldy- Wouthuysen 



and Pauli reductions of a Dirac hamiltonian containing a generic potential 
with harmonic time dependence. They found that differences do occur be- 
yond first order in 1/M. Detailed calculations show that these differences 
between the Pauli and Foldy- Wouthuysen reductions are small when the nu- 
clear potentials are ignored in the nuclear current operators. This seems 
the only consistent way to compare the operators since we use two different 
hamiltonians for the Pauli calculations. 

We discuss below the effects that the presence of the potentials in the 
non-relativistic current operators have on calculated observables. 



6 Results of Non-Relativistic Calculations 

We now discuss results of numerical calculations using the proper non-rela- 
tivistic reduction presented above. In Fig. 2. we show results for the reaction 
on an ^^O target with the same kinematics as in Fig. 1. Figure 2(a) shows 
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the spectral function while 2(b) shows the proton polarization. The non- 
relativistic calculations show the same effects due to the inclusion of the nu- 
clear potentials in the interaction operators that we saw in the corresponding 
convergence calculations of Fig. 1. The first and second order calculations 
without potentials in the nuclear currents ( dotted and dot-dashed curves 
respectively ) yield very similar results. This is generally true in the cases we 
have considered; going from first to second order in 1/M does little to move 
the results in the direction of the relativistic calculations. When potentials 
are included in the nuclear currents, a large change is noticeable in going 
from first order to second order calculations, particularly at larger values of 
missing momentum. Note that the non-relativistic calculations for the spec- 
tral function converge to a lower value than the simple expansion in powers of 
(E + M)~^ ( i.e. below the relativistic calculation ). This is because the nor- 
malization of the Schrodinger-cquivalcnt bound state wave function to unity. 
This results in the non-relativistic expansion converging at a point which is 
not the relativistic one, but a factor of the square of the inverse normalization 
constant lower than the fully relativistic result. This amounts to a reduction 
of the spectral function from the relativistic result by a factor typically in 
the range 1.2 to 1.4. Spin observables are not affected by changes in overall 
normalization, so the proton polarization calculations shown in Fig. 2(b) are 
very similar to those shown in Fig. 1(b), with slight differences coming from 
the replacement E ^ M in the non-relativistic nuclear current operators. 

Figure 3 emphasizes the behavior of the spectral function for the high 
missing momentum region of Fig. 2(a), with the missing momentum in the 
range 150 MeV/c to 300 MeV/c. In this region the first and second order 
calculations without potentials lie above the relativistic calculations, while 
the inclusion of potentials in first and second order moves the results to 
he below the relativistic results. Note that the the relativistic calculations 
were fitted to the data in the low missing momentum region, but still do 
rather well for high missing momenta. Similar results are obtained for "^^Ca 
and ^^Zr targets, whether the residual nucleus is left in the ground or ex- 
cited state. When the potentials are not included in the nuclear currents the 
results diverge from the relativistic calculations as the magnitude of the miss- 
ing momentum is increased. On the other hand, calculations which include 
potentials in the nuclear current operators remain close to the relativistic re- 
sults over a wide range of missing momenta. Note that we are only including 
terms to second order in the inverse mass. 
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Figure 4 shows non-relativistic calculations of the spectral function and 
proton polarization for the same reaction discussed in the previous figures, 
however in this case the energy of the incident electron is 2000 MeV, and 
the kinetic energy of the detected proton is fixed at 400 MeV. The larger 
energies allow for a much larger range of missing momenta than considered 
previously. It is important to note that the first and second order calcu- 
lations of the spectral function, without potentials included in the nuclear 
current operators, differ from the relativistic calculations by up to an or- 
der of magnitude for large missing momenta, while inclusion of the nuclear 
potentials results in convergence to the fully relativistic results in the high 
missing momentum region. In addition we see that for low missing momenta 
the convergence point is lower than the relativistic ( see insert ). The ( 7,p ) 
reaction shows behavior consistent with these observations for ( e, e'p ) at 
high missing momentum |TT|. The momentum transfer in the ( 7,p ) reaction 
is generally in the range 400 MeV/c to 600 MeV/c so these two reactions can 
both probe this part of the single particle bound state wave function. 

Proton polarization is shown in Fig. 4(b). In the region of large miss- 
ing momentum there are large differences between the polarization calculated 
with and without nuclear potentials in the current operator. The polarization 
calculated with first and second order currents containing nuclear potentials 
yields results close in magnitude and shape to the results of the fully rela- 
tivistic calculations. Note in particular, that in the region of the minimum 
and maximum in the relativistic calculations close to pm = — 400MeV"/ c and 
Pm = 4:00MeV/ c respectively, the calculations without potentials included do 
not reproduce the shape of the relativistic calculations at all. The potentials 
must be included in the nuclear current operators in order to get close to the 
relativistic results. In particular a measurement of the proton polarization 
near = — 400MeV^/c provides a clear opportunity to differentiate between 
relativistic and non-relativistic models. 

Calculation of the asymmetry parameter of Eq. (|T3p, for small values of 
angular momentum L of the bound nucleon ( L < 2 ), yields similar results for 
all the calculations whether fully relativistic; or non-relativistic first or second 
order, with or without potentials in the nuclear current operator. When the 
angular momentum of the bound nucleon is increased, the differences between 
these calculations of the asymmetry become larger, as is evident in Fig. 5. 
The asymmetry is calculated for a ^^Zr target, with the residual state in ^^Y 
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assumed to be a I/5 proton hole. The incident electron has an energy of 
461 MeV, and the kinetic energy of the detected proton is fixed at 100 MeV. 
In this case the differences are particularly apparent for missing momenta in 
the neighborhood of 20 MeV/c. 

7 Conclusions 

In order to clarify the differences arising from relativistic and non-relativistic 
descriptions of quasifree electron scattering [|, |[, we have discussed an ex- 
pansion of the S-matrix for the reaction ( e, e'p ) in powers of {E + M)~^ 
through the effective Pauli scheme. The resulting S-matrix depends on 
Schrodinger-like wave functions for the bound and continuum nucleons, and 
nuclear current operators which contain the strong Dirac potentials at the 
different orders. When the Dirac potentials are included in the nuclear cur- 
rent operators, the series essentially converges to the fully relativistic results 
at second order for light- to medium-mass nuclei we have considered. This 
indicates the importance of the role played by the nuclear potentials in the 
modification of the currents. When the potentials are not included in the 
nuclear currents, the calculations can be far from the relativistic results par- 
ticularly for larger missing momenta. 

These points were further studied in setting up a comparison between rela- 
tivistic and non-relativistic calculations. A proper non-relativistic calculation 
is obtained through several steps: normalization of the bound Schrodinger- 
like wave function to unity, proper normalization of the continuum Schrodinger- 
like wave function, and in the nuclear current operators the energy of the 
nucleons is set equal to the nucleon mass ( i.e. take the limit E ^ M ). 
An additional step of removing the Dirac potentials from the resulting nu- 
clear current operators yields the standard non-relativistic amplitude. This 
results in a consistent and fair comparison between the relativistic and non- 
relativistic calculations. The potentials used for the bound and continuum 
protons yield both the relativistic and non-relativistic wave functions, with 
normalizations handled appropriately. 

The non-relativistic calculations we have shown for first and second order 
nuclear current operators without potentials give the same results that a stan- 
dard non-relativistic calculation would give if provided with the Schrodinger 
equivalent wave functions derived from the Dirac equation. Inclusion of the 
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nuclear potentials in the non-relativistic nuclear current operator results in 
a large change in the calculated observables. In particular, calculations of 
the spectral function and final proton polarization using second order nu- 
clear current operators which include the Dirac potentials, can reproduce 
the magnitude and shape of the fully relativistic calculations. This is true 
even at large missing momenta where the non-relativistic calculations with- 
out potentials in the nuclear current operators yield very different results 
than the fully relativistic calculations. The polarization of the final proton is 
particularly sensitive to differences in the calculations, and measurements of 
this observable at large missing momenta could assist in the choice between 
the relativistic and non-relativistic approaches. 

We have also calculated the asymmetry defined in the text for the differ- 
ent orders, with and without potentials, and found that in cases in which the 
angular momentum of the bound nucleon is less than 2, there are no notice- 
able differences between these calculations and the full relativistic. When the 
orbital angular momentum of the bound nucleon is greater than 2 differences 
between the resulting asymmetry in missing momentum of these calculations 
will appear. This observable thus will be useful in differentiating between rel- 
ativistic and non-relativistic models only for nuclear states with large orbital 
angular momentum. 

Other groups have examined the sensitivity of the models to changes 
in the optical potentials and modifications of the wave functions ||, |], 
and have found sensitivities at the level of 15%. However, the essential 
differences between the relativistic and non-relativistic approaches do not lie 
in modifications of the wave functions. The essential difference comes from 
the appearance of the nuclear potentials in the nuclear current operators; 
a result of the reduction of the relativistic amplitude. We emphasize that 
these nuclear medium effects, characteristic of the present model, will not 
appear through a non-relativistic impulse description of the process. They 
are, however, inherent in the relativistic description. 
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Figure Captions 



FIG. 1. Observables for the reaction ^^0{ e, e'p where ^^N is in the Ipi 
state. The energy of the incident electron is 456 MeV, and the kinetic energy 
of the detected proton is fixed at 90 MeV with parallel kinematics. Hartree 



bound state wave functions are used [0 and the proton optical potentials 



are from [|I^ . The data are from reference . (a) spectral function and (b) 
proton polarization. Curves labelled according to their order in {E + M)~^ 
and whether or not the Dirac potentials are included in the nuclear current 
operators: dotted curve - first order in {E + M)~^ without Dirac potentials; 
dashed curve - first order with potentials; dot-dashed curve - second order 
without potentials; dot-dot-dashed curve - second order with potentials; solid 
curve - fully relativistic calculation. 

FIG. 2. Observables for the reaction ^^0{ e, e'p )^^A^ where ^^A^ is in a Ipi 
proton hole state. The kinematics are those of Fig. 1. Curves labelled ac- 
cording to their order in 1/M and whether or not the Dirac potentials are 
included in the nuclear current operators: dotted curve - non-relativistic cal- 
culations, first order in {E + M)~^ without Dirac potentials; dashed curve - 
first order with potentials; dot-dashed curve - second order without poten- 
tials; dot-dot-dashed curve - second order with potentials; solid curve - fully 
relativistic calculation. Potentials and data from the sources of Fig. 1. 

FIG. 3. Spectral function for ^^0{ e, e'p where ^^A^ is in the Ipi state. 
The kinematics are those of Fig. 1. Curves labelled as in Fig. 2. Potentials 
and data from the sources of Fig. 1. 

FIG. 4. Observables for ^^0( e, e'p Y^N where ^^N is in the Ipi state. The 
energy of the incident electron is 2000 MeV, and the kinetic energy of the 
detected proton is fixed at 400 MeV with parallel kinematics. Curves labelled 
as in Fig. 2. Potentials from the source of Fig. 1. 

FIG. 5. Asymmetry in missing momentum for ^^Zr{ e, e'p )^^Y where ^^Y 
is in the l/s state. The energy of the incident electron is 461 MeV, and 
the kinetic energy of the detected proton is fixed at 100 MeV with parallel 
kinematics. Curves labelled as in Fig. 2. Potentials from the source of Fig. 
1. 
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Figure 1 
M. Hedayati-Poor 



22 




-300 -200 -100 100 200 

missing momentum ( MeV/c ) 



Figure 2 
M. Hedayati-Poor 



23 




160 180 200 220 240 260 280 300 
missing momentum ( MeV/c ) 



Figure 3 
M. Hedayati-Poor 



24 




-600 -400 -200 200 400 

missing momentum ( MeV/c ) 



Figure 4 
M. Hedayati-Poor 



25 




Figure 5 
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